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Spatial clustering of USS sources and galaxies*. 
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ABSTRACT 

We present measurements of the clustering properties of galaxies in the field of redshift 
range 0.5 < z < 1.5 Ultra Steep Spectrum (USS) radio sources selected from Sydney 
University Molonglo Sky Survey and NRAO VLA Sky Survey. Galaxies in these USS 
fields were identified in deep near-IR observations, complete down to K s — 20, using 
IRIS2 instrument at the AAT telescope. We used the redshift distribution of K s < 20 
galaxies taken from Cimatti et al. (2002) to constrain the correlation length rrj. We 
find a strong correlation signal of galaxies with K s < 20 around our USS sample. A 
comoving correlation length ro = 14.0 ± 2.8 hr 1 Mpc and 7 = 1.98 ± 0.15 are derived 
in a flat cosmological model Universe. 

We compare our findings with those obtained in a cosmological N-body simu- 
lation populated with GALFORM semi-analytic galaxies. We find that clusters of 
galaxies with masses in the range M = io 13 - 4 ~ 14 - 2 hr 1 Mq have a cluster-galaxy 
cross-correlation amplitude comparable to those found between USS hosts and galax- 
ies. These results suggest that distant radio galaxies are excellent tracers of galaxy 
overdensities and pinpoint the progenitors of present day rich clusters of galaxies. 
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1 INTRODUCTION 

In hierarchical galaxy formation models, cosmic structures 
form by the gravitational amplification of small primor- 
dial fluctuations of matter density in the early Universe 
JWhite fc Reeslll97SD . Studies of the clustering properties 
of galaxies at high redshift are essential for understanding 
galaxy and structure formation. One of the most commonly 
used statistics to measure the clustering of a population 
of sources is the two-point correlation function which 
measures the excess probability of finding a pair of objects 
at a separation r with respect to a random distribution. 

High redshift radio galaxies are ideal targets for pin- 
pointing massive systems. Radio galaxies follow a close 
relation in the Hubble K — z diagram llLillv fc Longairl 
1984). The nature of this behavior in the K — z diagram 



* Based on observations obtained with the Australia Telescope 
Compact Array, the Anglo- Australian Telescope (Program 70. A- 
0514) 



shows that the stellar luminosities of 2 > 1 radio galax- 
ies are more luminous than normal galaxies at these red- 
shifts jDe Breuck et all 120021. At lower redshi fts they are 
found frequently in m oderately rich clusters jHill fc Lilly! 
ll99ltlYates et allll989f) Recently, galaxy overdensities com- 
parable to that expected for clusters of Abell class rich- 
ness are found near radio galaxies up to z = 1.6 feestl 
1200(1 iBest et all 12001 iBornancini et al.ll2004h . There is in- 
creasing evidence that galaxy overdensities around radio 
galaxies existe d at very hi g h red s hifts. Using Lyq and /o r 
Ha techniques. | Kurk et all i2000T ): fVenemans et al] (2002); 
iMilev et al] ll2004f) found overdensities of companion galax- 
ies around powerful 2 < z < 4 radio galaxies. 

In this paper we estimated the spatial correlation length 
for galaxies in the fields of USS targets selected from the 843 
MHz Sydney University Molonglo sky Survey (SUMSS) and 
1.4 GHz NRAO VLA Sky Survey (NVSS) in the redshift 
range 0.5 < z < 1.5, through the Limber's equation using 
an appropriate observed redshift distribution. We compared 
our results with those obtained in cosmological N-body sim- 
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ulations. This paper is organized as follows: Section 2 de- 
scribes the sample analyzed, we investigate the USS-galaxy 
cross-correlation analysis in Section 3. In Section 4 we in- 
terpret our results with those obtained in a cosmological N- 
body simulation. Finally we discuss our results in Section 
5. 

Throughout this paper we will use a flat cosmology with 
density parameters Qm = 0.25, Qa = 0.75 and a Hubble 
constant Ho = 100 h km s" 1 Mpc -1 . 



2 THE DATA 

The USS sample selection, radio data and redshifts 
used for this a n alysis was p resented and described by 
iDe Breuck et all (|2004l . 120051) . Detailed descriptions of 
the construction of th e galaxy catalogue is given in 
iBornancini et ail J2005). In summary, we used 11 K s - 
band images centered in Ultra Steep Spectrum (USS) ra- 
dio sources selected from the Sydney University Molonglo 
sky Survey and NRAO VLA Sky Survey, obtained with the 
instrument IRIS2 at the AAT telescope. We selected the 
redshift of the USS sample in the range 0.5 < z < 1.5. The 
redshift upper limit was adopted in order to assure that 
Tf-band images are sufficiently deep while the lower limit 
precludes the use of too close USS for which the field of 
view subtends too small linear scales. 

Our sample sources is listed in Table 1 in IAU J2000 
format, including the K s counterpart magn itude and the 
spectroscopic redshift iDe Breuck et alj 2005). 



3 USS-GALAXY CROSS-CORRELATION 
ANALYSIS 

The spatial USS-galaxy cross-correlation function £ us (r) is 
defined as the excess probability dP of finding a galaxy in 
the volume element dV at a distance r from a USS target, 



dP = ft [1 + &,,(!•)] dV, 



(i) 



where n is the mean space density of galaxies. The spatial 
two-point cross-correlation function £ U g{r) has been shown 
to be well approximated by a power law: 

In order to obtain the cross-correlation length ro 
we first determine the projected cross-correlation function 
Wug(er), where a is the projected separation between a USS 
target and a galaxy at redshift z. 

We use the follo wing estimato r of the projected cross- 
correlation function, <Peebledll98ch : 



n R UG(a) 



1, 



(3) 



n G DR(a) 

where tig and ur are the numbers of galaxies in the sample 
and in a random sample respectively, UG(a) is the number 
of real USS-galaxy pairs separated by a projected distance 
in the range a, a + 8a, and DR(a) are the corresponding 
pairs when considering the random galaxy sample. We esti- 
mate the corresponding correlation length using the Limber 
equation iLimberlll95ot) . The power law model for £u S (r) 
gives: 
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Figure 1. Distance distribution for K s < 20 galaxies taken the 
K20 survey iCimatti et alj |2002 ) . The dashed curve shows the 
selection function obtained fitting a Gaussian function. The upper 
small box shows the corresponding redshift distribution of the 
sample. 



(a) = By/lT 



r[(7-i)/2] 



r(7/2) 



(4) 



where the constant B regulates the amplitude of the cor- 
relation function taking into account the differences in the 
selection func tion of USS targets and galaxies and can be 
calculated by IlLilie fc Efstathioul ll988'): 



B 



(5) 



where N(jji) is the selection function of the galaxy survey, yi 
is the distance to USS target i and the sum extends over all 
USS targets in the sample. Equation 0] can be easily solved 
analytically if we perform a linear interpolation of uo ug be- 
tween its values at the measured er's. In order to calculate 
the constant B we evaluated the selection function N(yi) 
using the redshift (spectroscopic and photometric) distrib u- 
tion of K s < 20 galaxies published bv lCimatti et all ((2002) 1 . 
We model the distance distribution using a Gaussian func- 
tion, and find that y = 2000 h^ 1 Mpc (mean) and a y = 635 
h^ 1 Mpc (standard deviation) is a reasonable set of values 
that reproduces the distribution (See Figure . In Figure 
[5] we show the projected cross-correlation function w ttg (cr) 
for USS targets with spectroscopic redshifts in the range 
0-6 < z < 1.5 and galaxies with K s < 20. We estimate cross- 
correla tion functio n error bars using the jackknife tech- 
nique jEfrorJll982l) . We find a comoving correlation length 
r = 14.0 ± 2.8 /i" 1 Mpc with slope 7 = 1.98±0.04. 

We tested the accuracy of this result by varying the 
N(yi) distribution over a reasonable range of values for y and 



1 Data and further information available 
http : //www . arcetri . astro . it/~k20/releases 



at 
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Figure 2. Projected cross-correlation function u} ug (o) for USS 
targets with spectroscopic redshift in the range 0.6 < z < 1.5. 
The error bars correspond to the 1 a uncertainty estimated using 
the jackknife technique. 

a y . We find that the calculated ro is only weakly dependent 
on the y used, and is only affected at the 10% level. 



4 COMPARISON WITH N-BODY 
SIMULATIONS 

We interpret our results with the aid of a cosmological N- 
body simulatio n populate d with GALFORM semi-analytic 
galaxies JCole et alJl200rJl at different outputs correspond- 
ing to different redshifts, z = 0,1, and 3. This simulation 
was kindly provided by the Durham group. The cosmologi- 
cal model corresponds to matter and cosmological constant 
density parameters Q m ~ 0.25, S1a = 0.75, a power spec- 
trum tilt n = 0.95, an amplitude of fluctuations of as = 0.8, 
and a Hubble constant of Ho = 100 h kms~ Mpc -1 , where 
h = 0.7. The total number of particles is 1080 3 , the mass res- 
olution is 5.05 x 10 10 1i -1 Mq, and the number of dark matter 
haloes with masses greater than M = 10 12 h _1 MQ ranges 
from ~ 400, 000 at z = 3 to ~ 2, 200, 000 at z = 0. The 
number of GALFORM galaxies ranges from ~ 20, 000, 000 
to 120, 000, 000 at z = 3 and 2 = respectively. 

We now briefly explain the procedure by which galax- 
ies are assigned their properties in the Semi-analytic code. 
GALFORM is run for each halo in the numerical simula- 
tion, where galaxies are assigned to a randomly selected 
dark matter halo particle. The different galaxy properties 
such as magnitudes in different bands, including the K— 
band, depend on the dark matter halo merger tree. This 
merger tree is generated via Monte-Carlo modelling based 
on the extended Press-Schechter theory, and the evolution 
of the galaxy population in the halo is followed through time 
and different processes are considered in this evolution, in- 
cluding gas cooling, quiescent star formation and star forma- 
tion bursts, mergers, galactic winds and super winds, metal 



enrichment, extinction by dus t. For full d e tails o n the mod- 
elling the reader is referred to lCole et alJ <l2000lL 

We calculate the cross-correlation function using the 
simulation haloes with masses above a lower mass limit as 
centres, and as tracers, the GALFORM semi-analytic galax- 
ies. By comparing these measurements to the results from 
the cross-correlation between USS and normal galaxies, we 
make the implicit assumption that USS galaxies reside at the 
centres of dark-matter haloes. This comparison will make it 
possible to infer the mass of the structures associated to the 
USS hosts. 

Figure[3]shows the resulting real-space cross-correlation 
functions between haloes and semi-analytic galaxies at z = 1 
(top panel) for different halo masses. The shaded area cor- 
responds to the power law fit for the real-space correla- 
tion function inferred from the cross-correlation function 
for USS radio sources with spectroscopic redshifts in the 
range 0.6 < z < 1.5 (See Figure 2). In the middle panel 
of this figure, we compare the values of USS-galaxy cross- 
correlation length as a function of halo mass for three dif- 
ferent redshift outputs from the numerical simulations; as 
can be seen, the observed values are consistent with cluster 
masses within M = lo 13 4 " 14 2 h' 1 M at redshift z=l, in- 
dicating that our USS sample resides in massive clusters. In 
order to check whether our observational estimate of ro is 
affected by systematic biases, we calculate the projected cor- 
relation function in the numerical simulation and recover the 
real-space correlation length using Eq.^] setting B — 1; we 
consider the same range of separations available in the real 
data. The results for z — 1 are shown in filled circles in the 
middle panel. As can be seen, our conclusions on the mass 
of USS host haloes changes only slightly to M = io 13 - 2 " 13 - 8 
h^ 1 Mq, although we note that for lower and higher halo 
masses, the value of ro recovered from projected correlations 
is underestimated and overestimated, respectively. This test 
provides a useful check of our observational results which 
were derived from relatively small projected scales (r p <1 
Mpc). Our findings in the simulations indicate that reliable 
ro values are obtained using the power law approximation 
applied to projected correlations for r p < l/i -1 Mpc when 
the true correlation length is lesser than ~ 15/i _1 Mpc cor- 
responding to host halo mass M < 10 14 Mq. 

A further indication of the mass of USS galaxy host 
haloes comes from the lower panel of this figure, where 
the lines correspond to the projected cross-correlation func- 
tion measured in the GALFORM simulation for different 
masses (High to low masses from top to bottom lines at 
log 10 (cr/h _1 Mpc) = —0.3). H(cr) is calculated directly us- 
ing, 



S(ff) = Norm. 2 / £( vV 2 + 7r 2 )d7r, (6) 
Jo 



where we have used ir ma x — 80 h Mpc, and the nor- 
malisation, Norm., is set so that E(cr) and w(c) coincide at 
Iog 10 (o-/h _1 Mpc) = -1. The gray area shows the measured 
values of w(er) from the USS sample; as can be seen the mea- 
sured projected correlation function is in best agreement for 
M ~ 10 13 - 85 h _1 M ffl . 
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Figure 3. Top panel: cluster-galaxy cross-correlation functions 
from the simulation output at Z. = 1 for increasing halo masses 
Log(M) = 12.35,12.68,13.2,13.85 and 14.57 /i" 1 M Q (bottom 
to top). The shaded area shows the real— space correlation func- 
tion measured from the USS targets with spectroscopic redshift 
in the range 0.6 < z < 1.5. Middle panel: the real-space cor- 
relation length, ro, as a function of halo mass for three differ- 
ent redshift outputs; the open symbols show the values recovered 
from the real-space correlation function, whereas the filled sym- 
bols (shown only for 2 = 1) indicate the values of ro obtained 
from the projected correlation function following a similar proce- 
dure to that applied to the real data. The shaded area shows the 
allowed range of ro for the USS results. Bottom panel: the pro- 
jected USS cross-correlation function (gray area) and 3(<j) from 
the GALFORM simulation for the same halo masses as in the top 
panel; the normalisation is set so that 2(<r) and oj(ct) coincide at 
log 10 (<7/h- 1 Mpc) = -1. 



galaxies. Our testing with the numerical simulations also 
indicate that these observational results are not severely 
affected by the relatively small projected scales explored. 
We notice that for larger spatial correlation lengths, the 
power-law extrapolation for observed projected correlations 
in ~ l/i _1 Mpc fields would not give confident results. 

Previous studies have addressed the clust ering of galax- 
ies around radio sources. IWold et alJ ( 200o|) , using differ- 
ent cosmological parameters obtained a radio quasar-galaxy 
cross-correlation length about twice as large as the local 
galaxy autocorrelation length in suit able agreem e nt wit h our 
findings. The more recent work by iBarr et ail |2003) also 
indicates a moderate galaxy density enhancement around 
radio loud quasars similar to Abell richness clusters. 

Our analysis suggest that distant luminous radio galax- 
ies are excellent tracers of galaxy overdensities and may pin- 
point the progenitors of present day rich and moderate rich 
clusters of galaxies. 
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5 CONCLUSIONS 

We have analyzed the clustering properties of galaxies in 
the field of 0.5 < z < 1.5 Ultra Steep Spectrum (USS) ra- 
dio galaxies selected from the Sydney University Molonglo 
Sky Survey (SUMSS) and NRAO VLA Sky Survey (NVSS). 
We estimated the spatial clustering correlation length for 
galaxies in these fields, using the Limber equation using an 
appropriate observed redshift distribution, and we examined 
the dependence of galaxy clustering on USS targets luminos- 
ity. A comoving correlation length r = 14.0 ± 2.8 h~ x Mpc 
is derived and a slope 7 = 1.98±0.15. 

From our comparison with numerical simulations, we 
find that clusters of galaxies with masses in the range 
M = 10 13 ' 4 - 14 2 f l ~ 1 Mq have a cluster-galaxy correlation 
amplitude comparable to that found between USS hosts and 
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Table 1. USS Sample characteristics. Designation in IAU J2000 format, together with the K s counterpart magnitude, and the 
spectroscopic redshift (De Breuck et al, in prep.). 



(1) 




(2) 


(3) 


Name 




K— mag. 


z 






MAG_BEST 


spectroscopic 


NTVGJQ Tm ^9^9 
IN V DO JUlOZOZ- 


QQQ0^9 


i a 9fi_i_n ri9 


U.D14o±U.UUl 


IN V Do JU10044- 


"OOUDOO 


ID.yoztU.UO 


l.U4to±U.UUz 


NVSS J021716- 


-325121 


18.75±0.20 


1.384±0.002 


NVSS J030639- 


-330432 


17.88±0.11 


1.201±0.001 


NVSS J202026- 


-372823 


18.56±0.15 


1.431±0.001 


NVSS J204147- 


-331731 


16.86±0.05 


0.871±0.001 


NVSS J225719- 


-343954 


16.53±0.02 


0.726±0.001 


NVSS J230203- 


-340932 


17.34±0.07 


1.159±0.001 


NVSS J231519- 


-342710 


18.10±0.13 


0.970±0.001 


NVSS J234145- 


-350624 


15.94±0.04 


0.641±0.001 


NVSS J234904- 


-362451 


17.63±0.14 


1.520±0.003 
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